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Three-body systems with large scattering length display universal phenomena associated with a 
discrete scaling symmetry. These phenomena include resonant enhancement of three-body loss rates 
when an Efimov three-body resonance is at the scattering threshold. In particular, there can be 
resonant peaks in the atom-dimer relaxation rate for large positive scattering length. We improve 
upon earlier studies and calculate the atom-dimer relaxation rate as a function of temperature 
using a Bose-Einstein distribution for the thermal average. As input, we use calculations of the 
atom-dimer scattering phase shifts from effective field theory. 
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Introduction: The Efimov effect is a phenomenon that 
occurs in a nonrelativistic 3-body system with resonant 
short-range interactions. Efimov discovered in 1970 that 
there are infinitely many 3-body bound states with an 
accumulation point at the scattering threshold when the 
S-wave scattering length a is tuned to the unitary limit 
l/a = 0: 



(n) 



(1) 



where m is the mass of the particles, sq ~ 1.00624, and 
the binding wavcnumber of the branch of Efimov 
states labeled by n,. The geometric spectrum in ([T)) is 
the signature of a discrete scaling symmetry with scaling 
factor e'^/'*" ~ 22.7. For a finite scattering length that 
is large compared to the range of the interaction, the 
universal properties persist but there will only be a finite 
number of Efimov states. The corrections to the unitary 
limit are calculable in perturbation theory 0, 0, 0] ■ The 
Efimov effect is just one example of universal phenomena 
in the 3-body system with large scattering length 
For reviews of this "Efimov physics" , sec Rcfs. 
Here we focus on identical bosons in a single spin state 
and positive scattering length. 

Since we are interested in applications to ultracold 
atoms, we refer to the bosons as atoms, their 2-body 
bound states as dimcrs, and their 3-body bound states 
as trimers. If the scattering length a is large and posi- 
tive, there is a shallow dimer with binding energy = 
h"^ /{ma?). Moreover, there is an infinite sequence of pos- 
itive values of a for which there is an Efimov trimer at 
the atom-dimer scattering threshold HI: a = (e'^/'*°)"aH., 
where a, ~ 0.0708k~^ Another example of Efimov 
physics is an infinite sequence of positive values of a for 
which the 3-body recombination rate into the shallow 
dimer vanishes [ill, [13, [i3l ■ The universal aspects of 
Efimov physics are determined by two parameters: the 
scattering length a and the Efimov parameter k*. 

The alkali atoms used in most ultracold atom experi- 
ments have many deeply-bound diatomic molecules (deep 
dimers). Efimov physics is modified by the existence of 
the deep dimcrs because they provide inelastic channels 



for scattering processes and the decay of trimers. If there 
is an Efimov trimer near the scattering threshold, there 
will be a resonant enhancement of inelastic scattering 
processes. When a w (e'^/*'')"a*, the resonant inelastic 
process is dimer relaxation, in which the collision of an 
atom and a shallow dimer produces an atom and a deep 
dimer If there arc deep dimcrs, the universal phe- 

nomena associated with Efimov physics are determined 
by three parameters: a, k*, and a parameter 77* that de- 
termines the widths of Efimov trimers [l^ . 

Experimental evidence for an Efimov trimer in ultra- 
cold Cs atoms was recently provided by their signature 
in three-body recombination rates [l^. This signature 
could be unravelled by varying the scattering length a 
over several orders of magnitude using a Feshbach reso- 
nance. More recently, possible evidence for an Efimov 
resonance was also obtained in atom-dimer scattering 
[iHl, in three-body recombination in a balanced mixture 
of atoms in three different hyperfine states of ^Li [TtI . [l8l | , 
and in a heteronuclear system of "'^K and ^^Rb |19| . 

In this paper, we focus on the Efimov resonance in 
atom-dimer scattering found in [l6l |. This process was 
previously considered in [3, [13] ■ We go beyond these 
earlier studies in various respects: we use the full effec- 
tive field theory results for the atom-dimer phase shifts 
instead of the effective range expansion and perform a 
thermal average using the Bose-Einstein instead of the 
Boltzmann distribution. Moreover, we correct an error 
in the calculation of [20j . 

Atom-dimer scattering: We consider the scattering of 
an atom with mass = m and dimer with mass 
Too = 2m. The wavenumbers and po of the incom- 
ing atom and dimer, respectively, can be decomposed into 
the total wavcnumber ptot = Pa + P_d and the relative 
wavenumbcr k = fpA ^ \vd- Because of Galilei invari- 
ance, the scattering observables depend on the relative 
wavenumbers and the collision energy E in the center-of- 
mass system only: E = 3?i^fc^/(4TO) where k = |k|. The 
differential cross section for clastic atom-dimer scattering 
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where fADik,9) is the scattering amphtude. The elastic 
cross section cr^]^''*"^^ is obtained by integrating Eq. ([2]) 
over the full solid angle. The total cross section (includ- 
ing elastic and inelastic contributions) can be calculated 
using the optical theorem: 



(total) 
^AD 



An 



lmfADik,9 = 0), 



(3) 



such that the inelastic cross section is given by the dif- 
ference of the total and elastic cross section. At low en- 
ergies, higher partial waves with L > are suppressed 
and the scattering amplitude is dominated by S-waves 

{L = oy. 



fAoik) = [fccot(5ff^(fc)-zA:] 



(4) 



For the S-wave atom-dimcr phase shift kcotSQ^{k), we 
will use the results from a calculation using the effective 
field theory of Ref. [2l| . A convenient paramctrization of 
these results was given in 0, Q : 



ka cot S^^{k) 



where 



ci(fca) 
C2{ka) 
(j){ka) 



~ ci(fca) 

4- C2(fca) cot[so ln(0.19 a/o 

= -0.22 + 0.39fc2a2-0.17fc 

= 0.32-t-0.82fc2a^ - 0.14fcV, 

= 2.64- O.SSfc^a^ -h0.23fc ' 



(5) 
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This paramctrization is valid up to the dimer breakup 
wavenumber of fc^) = 2/(V3a). A paramctrization for 
higher wavenumbers beyond the dimer breakup exists 
[2 2( 1 but will not be required for our purposes as we will 
demonstrate below. 

To leading order in the large scattering length, atom- 
dimer relaxation can only proceed via S-waves. For the 
relaxation into deep dimers to take place, the atom and 
the dimer have to approach each other to very short dis- 
tances. However, because of the angular momentum bar- 
rier this can only happen in the relative S-wave. The 
parametrization of the S-wave phase shift in Eqs. (0 E]) 
is therefore sufficient to calculate atom dimer-relaxation. 

Atom-dimer relaxation: To incorporate the effects of 
deep dimers, we make the simple replacement in the am- 
plitude [13 



ln( 
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where 77* determines the probability for an atom and a 
dimer to scatter into an energetic atom and deep dimer 
pair at short distances. This inelastic process generates 
the width of the Efimov resonances. The phase shift be- 
comes imaginary even below the dimer breakup threshold 



and the released binding energy is converted to the ki- 
netic energy of the recoiling atom and dimer. They thus 
are lost to the system. The event rate f3 for this dimer re- 
laxation process in an ultracold gas of atoms and dimers 
can be written as 



dt 
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where and njj denote the number density of the atoms 
and dimers. respectively. 

For an ensemble of atoms and dimers at nonzero tem- 
perature that arc held in a trap, temperature and trap 
geometry have to be included in the calculation of the ob- 
served dimer losses. The dimer loss rate can be expressed 
as 



dh 



n 

i=A,D 



m, (9) 



where Nd is the number of dimers and we use the gen- 
eralized Bose-Einstein distribution function 



exp 
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with i = A, D denoting an atom or dimer, respectively. 
The properties of the trap enter via the average trap 
frequency uj while the function g(k) to be averaged is 
given by 
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In the limit fc ^ 0, g{k) reduces to the relaxation rate 
constant at zero temperature. Note that the function 
g{k) introduces an implicit dependence on the angle be- 
tween Pa and p/j. The chemical potentials fii are fixed 
via the equation 



d^r 



d^Pr 
(27r)= 



(12) 



with Ni being the particle number and i — A, D. 

All angular integrations except for one can be carried 
out immediately and the expression ^ can be rewritten 
as: 

^""^ = -2^1 ''"^'i P'otdptotJ^ k^dk 

X j dxnA{pA,r)nD{pD,r)g{k) , (13) 

where x is the cosine of the angle between ptot and k. 

In evaluating Eq. (fT^ . we will cut off the integral over 
k at the breakup wavenumber ku ^ 2/{^/ia) since the 
parametrization in Eqs. (OEl) is only valid up to k^. We 
have estimated the error from this simplification by using 
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FIG. 1: The dimer relaxation coefficient /9 as a function of 
a/ao for T — 170 nK, a, — 397ao, and different values of 
77,. The data points are from [l^. BO indicates a Boltzmann 
average. 



the unitary bound i/k for the S-wave scattering ampli- 
tude fAoik). For the data of Ref. the error involved 
is largest for the largest scattering length considered but 
even there only adds up to 0.2% for the largest temper- 
ature T = 170 nK. As a consequence, we can simply 
neglect the contribution from k > k]j in the analysis of 
the data. 

Results and Discussion: We now apply our formalism 
to the experimental data for the atom-dimer relaxation 
rate of ultracold Cs atoms as a function of the scattering 
length obtained by Knoop et al. [31 . Our free parame- 
ters are a, which determines the position of the resonance 
and rjif which determines its width. These parameters 
cannot be calculated in our approach and must be taken 
from experiment. We will determine a* and 77* from the 
data of Knoop et al. and compare our results with what 
is known from other experiments. In Ref. [iGj . the dimer 
relaxation coefficient was extracted from the dimer loss 
data using a loss model resulting in the rate equation 



dt 



-N, 



D 



PnAN£i+ dimer loss term , (14) 



with fiA = [mU]^ / {^TrkBTyf N A the mean atomic den- 
sity. In order to compare our calculation with the exper- 
iment of Knoop et al., we extract a value for (3 from our 
result for dNo/dt (cf. Eq. HSl)) using 



(3: 



27 d 



SuaNd dt 



D ■ 
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We start with the data at T = 170 nK and fix the 
chemical potentials as described above. For atom num- 
ber Na = 10^, dimer number Nd = 4 x 10'^, and an 
average trap frequency of tZ; = 45 Hz [23| . we obtain 
the chemical potentials ija = —2.74 x 10~^ fcsK and 
Hd = -8.17 X IQ-"^ ksK. In Fig. [1] the data for the re- 
combination constant /? is shown together with our best 



fit as the solid (blue) line. We only take into account 
data points for a > BOOao. This fit yields x^/dof = 1.2. 
We obtain for the peak position a» = 397ao and for the 
resonance width parameter 77* = 0.034. Also shown as a 
dashed-dotted (green) line is the resulting curve for the 
same resonance position but with = 0.06. This value 
of 7]^ was obtained from fitting the three-body recombina- 
tion resonance in Cs for negative scattering length in . 
It is also compatible with the three-body recombination 
data for positive scattering length presented in the same 
paper. ^ However, the data for a > are not very sensi- 
tive to the precise value of ry* and values of ry* as small 
as 0.01 would also be compatible. The width parame- 
ter 77* should only be weakly dependent on the magnetic 
field in a universal region |7[. In first approximation it 
can be assumed to remain constant. A more serious puz- 
zle is that the resonance position a* ~ 397ao extracted 
from the dimer relaxation data is not compatible with 
the value a* « 1200ao found in [l^. The three-body re- 
combination data cannot be satisfactorily described using 
a* = 397ao. This disagreement requires further study. 
However, one has to keep in mind that the atom-dimer 
resonance is at the border of the universal region since 
the van der Waals length scale is Ivdw ~ 200ao for Cs 
atoms. 

For temperatures much larger than the Bose-Einstein- 
condensation temperature, the thermal average can be 
replaced by a Boltzmann average. The critical temper- 
atures are estimated by setting the chemical potential 
to zero and solving Eq. for the critical temperature. 
This yields T^a « 94 nK for the atoms and T^d ~ 32 nK 
for the dimers. Therefore, it seems justified to use Boltz- 
mann distributions instead of the Bose-Einstein distribu- 
tion functions Ui . The resulting calculation for (3 is analo- 
gous to the method of Ref. [l^l but uses the parametriza- 
tion of the scattering phase ([5]) instead of an effective 
range expansion. The resulting curve is also shown in 
Fig. [B as the dashed (red) line. We remark that in the 
numerical evaluation of Eq. (8) in Ref. [l^l a factor of k^ 
in the thermal average was omitted. This lead to a wrong 
normalization of the curves in Figs. 1 and 2 of this refer- 
ence. The dashed line obviously also describes the data 
quite well but yields 77, = 0.036. This shows how tem- 
perature dependence and averaging methods can partly 
be compensated by adjusting the width parameter 77*. 
For an accurate determination of 77*, reliable tempera- 
ture and particle number measurements are thus crucial. 

Using the effective range expansion for the atom- 
dimer scattering amplitude as in Ref. [13 instead of the 
phase shift parametrization of Eqs. ^ ^ does not al- 
ter the overall shape or normalization of the dimer re- 
laxation coefficient /3. However, the extracted value of 



^ Note that in the Cs experiment of [TBI the regions a > and 
a < are not required to have the same parameters since they 
are separated by a zero in the scattering length rather than a 
pole. See also the discussion in [Tsll . 



4 




FIG. 2: The dimer relaxation coefficient /9 as a function of 
a/ao for T = 40 nK, a, = 397ao, and different values of 77,. 
The data points are from [l^ . 

a* is shifted by about 3% to a higher value whereas 
rjif remains unchanged. The scattering length approxi- 
mation with the atom-dimer scattering length given by 
o^adIo- = 1-46 + 2.15 cot[so ln(a/a») + ir\^\ does not lead 
to an equally good fit result. The obtained values for 
are smaller especially for higher values of a. Besides a 
change in peak position, width and height, we can only 
obtain a /dof of 3 — 4. 

We now turn to the data for T = 40 nK [l^ and 
compare them to our theoretical results for different 
values of 77*. Here, the atom and dimer numbers are 
Na = 3 X 10'* and Nd = 4 x 10'^, and the trap frequency 
is cU = 25.2 Hz [1^ leading to the chemical potentials 
^1A = -1-21 X 10-8 keK and fio = -8.87 x 10"* fesK. 
The critical temperatures are estimated as Tc.a ~ 35 nK 
for the atoms and Tc^d ~ 18 nK for the dimers such that 
the temperature is only slightly larger than the critical 
temperature for the atoms. In Fig. [2] we show the data for 
the relaxation coefficient /3 at T = 40 nK together with 
our fit results. The dashed (red) line gives our predic- 
tion for the relaxation coefficient /3 using the parameters 
obtained by fitting the 170 nK data. The prediction is 
about a factor 2 too large compared to the data. The 



dip in the data at the peak position represented by the 
third and fourth data point cannot be reproduced within 
our theory. If it is not simply a statistical fluctuation, it 
must be due to physics not captured in our theory such 
as non-universal effects or four-body physics [U, [2^, [2^ . 
If we keep the resonance position at a* = 397ao but fit 
the parameter 77* to the 40 nK data excluding the third 
and fourth data point, we obtain the solid (blue) line. 
Still excluding the third and fourth data point this gives 
77* = 0.016 with x^/dof = 1.5 and describes the data 
satisfactorily. 

In summary, we have calculated the atom-dimer relax- 
ation rate for large positive scattering length in a univer- 
sal zero-ranee approach. We have improved on previous 
studies [13, [2^ by using a Bose- Einstein distribution for 
the thermal average and calculations of the atom-dimer 
scattering phase shifts from effective field theory. Our 
results were then applied to the atom-dimer relaxation 
data for Cs obtained by Knoop et al. [31 . Fitting the 
resonance position and width parameters a* and 77*, we 
could get a good description of the 170 nK data. Us- 
ing these parameters, we overpredict the relaxation data 
at 40 nK by a factor of two. Moreover, our theory is 
not able to reproduce the dip in the 40 nK data directly 
at the resonance position and the corresponding physics 
appears to be missing in our theory. We demonstrated 
that this discrepancy is neither due to the thermal aver- 
aging procedure nor due to the phase shift parametriza- 
tion used. The resonance position at the border of the 
universal region and the mismatch in the extracted reso- 
nance position from atom-dimer relaxation and the three- 
body recombination data [Tsl suggests that non-universal 
physics could be responsible. However, it is also conceiv- 
able that four-body losses become important at the lower 
temperature. This question deserves further study. 
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